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Figure 1: The experimental setup of the NA49 experiment 



1. Introduction 



Elliptic flow has its origin in the spatial anisotropy of the initial reaction volume in non-central 
collisions and in particle rescatterings in the evolving system which convert the spatial anisotropy 
into a momentum anisotropy [[I]]. The spatial anisotropy decreases rapidly because of the fast 
expansion of the system making the momentum anisotropy measured at the end of this evolution 
strongly dependent on the matter properties and the effective equation of state (EoS) at the early 
stage [jsj, 0]. It is particularly sensitive to the degree of thermalization in the produced particle 
system. Flow of heavy particles is affected more strongly by changes in the EoS than flow of 
pions [Q, [Sj, 1^. Moreover, various hadron types are believed to decouple at different times and 
with different temperatures Thus the elhptic flow of various particle species allows insight into 
different stages of the collision. Comparison of measured anisotropics of different particle species 
with various model calculations, for example hydrodynamical or quark coalescence, provides an 
important test of various evolution scenarios. 

The anisotropic flow parameters measured to date at SPS and lower energies are mainly those 
of pions and protons [^, ^]. We have extended elliptic flow measurements in 158A GeV Pb+Pb 
(\/^NN = 17.2 GeV) collisions to A hyperons and particles to test the validity of the hydrody- 
namic scenario and check the degree of thermalization at SPS energies. 



2. Experiment and Data 

The NA49 experimental setup is shown in Fig. 1. The main components of the NA49 detector 



[|10p are four large-volume Time Projection Chambers for tracking and particle identification. The 
TPC system consists of two vertex chambers inside the spectrometer magnets which allow separa- 
tion of positively and negatively charged tracks and a precise measurement of the particle momenta 
with a resolution <j{p) jp^ = (0.3 — 7) x 10^"* (GeV/c)^^ Two main chambers, placed behind the 
magnets at both sides of the beam, were optimized for high precision detection of the ionization loss 
dEldx with a resolution of 3— 6%. Downstream of the TPCs a veto calorimeter detects projectile 
spectators and is used for triggering and centrality selection. The data sample consists of 3 x 10^ 
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Figure 2: The invariant mass distribution of pn^ candidate pairs in Pb+Pb collisions in the centrality range 
oIOtot = 19.5—23.5%. The polynomial parametrization of the background is indicated by the dashed hne. 
Red vertical lines represent invariant mass cuts. 



semi-central Pb+Pb events after online trigger selection of the 23.5% most central collisions. The 
events were divided into three different centrality bins, which correspond to the first three bins 
used in a previous analysis (see Table 1 in |^). They are defined by centrality ranges 0-5% (bin 
1), 5-12.5% (bin 2), and 12.5-23.5% (bin 3) which correspond to impact parameter ranges which 
are: 0-3.4 fm (bin 1), 3.4-5.3 fm (bin 2), and 5.3-7.4 fm (bin 3). The measurement in the centrality 
range alajoT = 5—23.5% (called mid-central) is obtained by averaging the results of bins 2 plus 3 
with weights corresponding to the fractions of the total cross section in these bins. 



3. Selection of particles 

The A hyperon candidates were selected from the sample of V°-track configurations consist- 
ing of oppositely charged particles, which include the A decays into proton and (branching 
ratio 63.9%). Geometric and quality criteria ensured that only reliably reconstructed tracks were 



processed. The identification method [11] relies on the evaluation of the invariant mass distribution 
and is enhanced by daughter particle identification applying a cut in dEldx around the expecta- 
tion value derived from a Bethe-Bloch parametrization. mesons, which decay into 71+ and 
(branching ratio 68.6%), were identified in a similar way as A hyperons. Instead of identification of 
daughter particles by dEldx measurement we assigned the proton mass to one of the daughter par- 
ticles and excluded those pair candidates which give entries near the expected A, A invariant mass 
peak, 1.100 < rripj^-, mpji+ < 1.132 GeV/c^ (compare Fig. 3 (left) and (right)). The yields of A 
hyperons and Kf particles are obtained by counting the number of entries in the invariant mass peak 
above the estimated background as a function of the azimuthal angle with respect to the event 
plane (see below). The background is estimated from a fit of sum of Lorentz distribution for or A 
and a polynomial background (Fig. 2 and Fig. 3). The number of accepted particles in mid-central 
events is about 740.00 for A hyperons in the invariant mass window 1.108-1.124 GeV/c^ and 
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Figure 3: The invariant mass distributions of n^n^ candidate pairs in Pb+Pb collisions at ^/snn = 17.2 GeV 
in the centrality range uIOtot = 19.5—23.5%. The distribution is with (left) and without (right) admixture 
of pTZ^, pn^ candidate pairs with an invariant mass near to the expected A, A hyperon peak (see text for 
details). The polynomial parametrization of the background is indicated by the dashed line. Red vertical 
lines represent invariant mass cuts. 
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Figure 4: Transverse momentum distributions for A hyperons (left) and particles (right) measured in the 
centrality range olojar = 19.5—23.5% of Pb+Pb collisions at ^/sNN = 17.2 GeV. 



440.000 for mesons in the invariant mass window 0.488—0.508 GeV/c^ . The acceptance of 
A hyperons covers the range 0.4 < /?t ^ 4.0 GeV/c (see Fig. 4(left)) and — 1-5 < 3^ < 1.0 and 
strongly depends on /jt and y. The range of transverse momentum of mesons is more narrow 
0-2 ^ /'T ^ 3.5 GeV/c (see Fig. 4(right)). Multiplicative factors were introduced for every A par- 
ticle to correct the A yields for detector and reconstruction efficiency which depends on pi and y. 
Such an additional correction has not been introduced for K'^ mesons so far. K*^ elliptic flow data 
are still very preliminary. 



4. Method 



The elliptic flow analysis is based on the standard procedure outlined in |12p to reconstruct 
the event plane for each event and the corrections for the event plane resolution. The event plane 
is an experimental estimator of the true reaction plane and is calculated from the azimuthal distri- 
bution of primary charged n mesons. Identification of pions is based on energy loss measurements 
(dE/dx) in the TPCs. To avoid possible auto-correlations, tracks associated with A or candidates 
are excluded from the event plane calculation. The method to determine the event plane angle ^2EP 
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Figure 5: Elliptic flow of A hyperons measured in mid-central Pb+Pb collisions as a function of rapidity. 
The open points have been reflected about midrapidity. 



uses the elliptic flow of pions, according to the formula: 



N 

^2 = £pV[co5(2Cb) - {cos{2(t>ub))], 

i=\ 
N 

Yl = £;?V[^'«(201ab) - (^i«(2^1ab))], 
i=l 



(4.1) 



<D2EP = tan-^ (^^) /2. 

where X2,Y2 are the components of the event plane flow vector Q2 and the sums run over accepted 
charged pion tracks. The acceptance correction is based on the recentering method of which 
consists of subtracting in Eq. 4.1 the mean values (co5(2^ab)) and (5/?i(2(^ab))- These mean values 
are calculated in bins of pj and rapidity for all charged pions in those events which contain at least 
one A hyperon or meson candidate. The means were stored in a 3-dimensional matrix of 20 pj 
intervals, 50 rapidity intervals, and eight centrality bins. 

A second level acceptance correction is done by using mixed events. We used 10 mixed events 
for each real event. Particles for mixed events are randomly selected from different events in the 
same centrality bin with at least one A or /T*' particle. The final angular distributions were obtained 
by dividing the real A and K'^ angular distributions by the mixed event distributions to remove the 
acceptance correlations remaining after recentering. The corrected particle azimuthal distributions 
are then fitted with a truncated Fourier series: 

dN 



fi?(01ab-<J>2Ep) 

const X (1 + vf'cos[2{^^-^2Ep)] 

+ vfW[4((/^ab-<J>2Ep)]). 



(4.2) 
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Figure 6: Elliptic flow of A hyperons (left) and mesons measured by NA49 in Pb+Pb and NA45/CERES 
in Pb+Au mid-central collisions as a function of transverse momentum. 
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Figure 7: Dependence of V2 on p-j for A hyperons (left) and mesons (right) measured at ^/snn = 17.2 
GeV and ^Jsnn = 200 GeV in comparison to a hydrodynamical calculation at ^Jsnn = 200 GeV for two 
centrality ranges. 



The elliptic flow V2 is evaluated by dividing the observed anisotropy v'^'^ by the event plane resolu- 
tion R: 



V2 



..ohs 

R 



The event plane resolution, 

/?=(c05[2(<I>2Ep-<I>r)]) 



2(co.[2(<I>^Ep-<I>'^Ep)]). 



(4.3) 



(4.4) 



is calculated from the correlation of two planes (*J*2EP'*^2Ep) random sub-events with equal 
multiplicity. The results are /? = 0.27, 0.34 and 0.40 for centrality bins 1, 2 and 3, respectively. 
The uncertainty due to the background subtraction, event plane resolution and the mixed event 
corrections are added to the total error. The observed hexadecupole anisotropy v"^^ is consistent 
with zero within statistical errors. 
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Figure 8: Plot of V2 of charged pions, protons, A hyperons and mesons as a function of pT in mid-central 
Pb+Pb collisions at ^Js^n =17.2 GeV. The solid lines represent results of a blast wave fit. 



5. Results 

Fig. 5 shows the /jt averaged elliptic flow obtained from all identified A hyperons without 
Pi cuts. It exhibits no significant dependence on rapidity as was also observed for protons (see 
Fig. 6 of [^) in mid-central events. The full rapidity range of the data in Fig. 5 was used to study 
V2 as a function of p^. The pj dependence of rapidity-averaged A and elliptic flow is shown in 
Fig. 6(left) and Fig. 6(right) respectively in comparison to CERES data[pj|]. The NA49 V2 param- 
eter significantly increases with transverse momentum and agrees with CERES results in Pb-i-Au 
mid-central collisions at the top SPS energy. Fig. 7 shows a comparison of V2{pt:) of A hyperons 
(left panel) and K*^ mesons (right panel) for mid-central events measured by the NA49 and STAR 



experiments []14[]. At SPS energy the elUptic flow grows Unearly with pj up to ~2 GeV/c, but the 
increase is more pronounced at RHIC than at SPS energy. It should be noted that RHIC mid-central 
data have been measured in the centrality range a/aTOT = 5—30% while SPS events are somewhat 
more central. The effect of different centrality ranges has been estimated by hydrodynamic calcu- 
lations [^] at RHIC energy for the slightly different centrality bins of NA49 and STAR. As shown 
by the corresponding curves in Fig. 7 this explains only partly the difference between both mea- 
surements. A comparison of V2(/?t) for pions, protons, A hyperons and mesons as measured by 
the NA49 experiment in mid-central events is displayed in Fig. 8. The values for pions and protons 
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Figure 9: Left: V2(pt) of charged pions, protons, A hyperons and K^^ mesons from 158A GeV Pb+Pb mid- 
central collisions. Hydrodynamic model calculations at ^/snn = 17.2 GeV are shown as solid (7y=120 MeV) 
and dashed {Tf=16Q MeV) lines (see text for details). Right: V2 scaled by the number of quark as a function 
of scaled pr- All data are from 158A GeV Pb+Pb mid-central collisions. The dashed line is the scaled result 
of the fit to p, A and K^. 



were obtained as the cross section weighted averages of the measurements published in for the 
appropriate centrality range. As seen in Fig. 8 the elliptic flow grows linearly with pj for all parti- 
cle species but the rise for pions starts from pj close to zero while for protons, A and mesons it 
starts from pj 0.5 GeV/c. The elliptic flow for pions is significantly larger than that for heavier 
particles although at ~ 2 GeV / c the flow becomes similar for all particle species. Except for 
mesons the measurements are reproduced by blast wave fits |Q, [l^, 17] (curves in Fig. 8) with 
the following parameters: freeze-out temperature T = 95 MeV, mean transverse expansion rapidity 
Po = 0.85, its second harmonic azimuthal modulation amplitude p2 = 0.021 and the variation in 
the azimuthal density of the source elements S2 = 0.035. In Fig. 9(left) the measured values of V2 
are compared to hydrodynamical model calculations [|l^] assuming a first-order phase transition to 
a QGP at the critical temperature Tc = 165 MeV. With the freeze-out temperature Tf = 120 MeV 
tuned to reproduce particle spectra, the model calculations significantly overestimate the SPS re- 
sults for semi-central collisions (full curves in Fig. 9(left)) in contradiction to predictions at RHIC 
energy which agree with data quite well for pj ^2 GeV/c [0]. The discrepancy at SPS may 
indicate a lack of complete thermalisation or a viscosity effect. On the other hand, the model 
reproduces qualitatively the characteristic hadron-mass ordering of elliptic flow. Thus the data sup- 
port the hypothesis of early development of collectivity. The calculation from the same model with 
a higher temperature Tf = 160 MeV exhibits better agreement with the A flow data (dotted curves 
in Fig. 9(left)). Unfortunately, the model does not simultaneously reproduce the mj spectra with 
such a high freeze-out temperature. A partial solution of this problem can be found by coupling a 
hadronic rescattering phase to the hydrodynamical evolution and hadronisation [Q]. 

Coalescence models predict that V2 will approximately scale with the number n of constituent 
quarks. When plotting at the similarly scaled value pj/n results for all hadrons are expected to 
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fall on an universal. This prediction agrees well with RHIC data in the intermediate transverse 
momentum region pj/n> 0.7 GeV/c [0]. The universal curve for all hadrons is expected to rep- 
resent the momentum-space anisotropy of constituent quarks prior to hadron formation. The naive 
coalescence model roughly agrees with our proton, A and K'^ at higher pj values ( Fig.9(right)) 
although the pj range of accepted particles is too narrow and error bars are to large to draw a clear 
conclusion. In more realistic quark coalescence model the resonance decays and quark momentum 
distribution in hadrons can lead to higher V2 for pions and generally deviations of meson elliptic 
flow from the scaling behavior as observed on Fig. 9(right) and also seen by the STAR experiment. 
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